Zero-energy states in corrugated bilayer graphene 
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Anomalous quantum Hall effects in single-layer and bilayer graphene are related with nontrivial 
topological properties of electron states (Berry phases it and 2tt, respectively). It was known that 
the Atiyah-Singer index theorem guarantees, for the case of the single-layer, existence of zero-energy 
states for the case of inhomogeneous magnetic fields assuming that the total flux is non-zero. This 
leads, in particular, to appearance of midgap states in corrugated graphene and topologically protects 
zero-energy Landau level in corrugated single-layer graphene. Here we apply this theorem to the 
case of bilayer graphene and prove the existence of zero-energy modes for this case. 

PACS numbers: 73.43. Cd, 02.40.Vh, 81.05.Uw 



I. INTRODUCTION 

Graphene, that is a two-dimensional allotrope of 
carbon formed by single carbon atom sheet is a subject 
of hot interest now (for review, see Refs. [l|0[3|[il) . 
One of the most interesting aspects of the graphene 
physics from theoretical point of view is a deep and 
fruitful relation with the quantum electrodynamics 
and quantum field t he or y 5 ->& 7 ->&£ i 1 1 1 1 1 1 2 . As was proven 
experimentally in Refs . [laflil charge carriers in the single- 
layer graphene are massless Dirac fermions characterized 
by "chirality", or "Berry phase" n. As a consequence, 
graphene demonstrates anomalous quantum Hall effect 
due to existence of zero-energy Landau level. The latter 
can be considered 2 - ' 11 ' 13 as a simple consequence of the 
famous Atiyah-Singer index theorem^ which plays an 
important role in the modern quantum field theory and 
theory of superstring a 16 ' 17 . 

Charge carriers in bilayer graphene, formed by 
two graphite atomic sheets 18 , can be, in a good 
approximation, considered as chiral fermions with the 
Berry phase 2tt which leads to another type of anomalous 
quantum Hall effec^^ and to a very unusual character 
of electron transmission through potential barriers 9 . 
Exact solution of the Schrodinger equation for the bilayer 
graphene in homogeneous magnetic field 1 ^ demonstrates 
existence of the zero-energy Landau level with twice 
larger degeneracy than for the case of single layer. 
However, topological origin of this feature was not 
clarified yet. This is the aim of the present work. We 
will prove that the existence of the zero-energy states 
in bilayer graphene is also a consequence of the Atiyah- 
Singer index theorem and thus is topologically protected. 

This is an important question since it is 
know n 20 ' 21 ' 22 ' 23 ' 24 ' 25 that graphene is always 
corrugated and covered by ripples which can be 
either intrinsi c 21 ' 22 ' 25 ' 26 or induced by a roughness of 
substrat o 23 ' 24 . In general, non-flatness of graphene leads 
to appearance of the pseudomagnetic inhomogeneous 
gauge fiel d 20 ' 27 acting on the charge carriers. Based on 
the topological arguments one can demonstrate that this 
pseudomagnetic field should also result in appearance of 



zero-energy states (pseudo-Landau levels) as was recently 
confirmed by model 28 and first-principles 2 ^, electronic 
structure calculations. This can provide a mechanism 
of formation of charge inhomogeneity in corrugated 
graphene 2 ^ and thus essentially effect on its electronic 
properties. Also, this "topological protectorate" of 
zero-energy Landau level can explain why it is narrower 
than the higher-energy levels as was recently observed 
experimentally 3 ^. The ripples in bilayer graphene has 
been already observed experimentally 2 - 2 * but it is still not 
clear whether this leads to formation of the zero-energy 
states, similar to the case of single layer, or not. Here we 
give a positive answer on this question based on a very 
general topological consideration. 



II. FORMULATION OF THE PROBLEM 

The corrugation leads to important consequences 
for the electronic structure of graphene. The nearest- 
neighbor hopping integral 7 turns out to be fluctuating 
due to its dependence on the deformation tensor— 
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where h is the displacement in the direction 

perpendicular to graphene plane, y) are 

coordinates in the plane and Ui are corresponding 
components of the displacement vector: 
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Taking into account this inhomogeneity in a standard 
tight-binding description of the electronic structure 
of graphene 4 - one can obtain an effective Dirac-like 
Hamiltonian describing electron states near the conical 
if-point: 



H = vpo (-ihV - A) 



(3) 



where vf = \/37oa/2fi and A is the "vector potential" 
connected with the deviations of the hopping parameters 



2 



7» from their unperturbed value 70: 

Al = ^—(72 + 73 - 27i). 
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(73 - 72) , 



(4) 



where the nearest neighbors with vectors 
(-o/V3, 0) ; (a/2y/H, -a/2) ; (a/2\/3, a/2) are labelled 
1,2, and 3, correspondingly, a is the lattice constant. 
This means that the flexural fluctuations act on the 
electronic structure near the if-point as an Abelian 
gauge field which is equivalent to the action of a random 
magnetic field. Thus, the bending of graphene violates 
the time-reversal symmetry for a given valley; of course, 
the Umklapp processes between K and K' points will 
restore this symmetry. As was suggested in Ref. [2(3 
these effective magnetic fields might be responsible for 
suppression of the weak localization effects in graphene. 

Whereas a smooth deformation of the graphene sheets 
produces the gauge field similar to electromagnetic one, 
different topological defects in graphene inducing inter- 
valley (Umklapp) processes can be considered as sources 
of a non- Abelian gauge field; corresponding analogy with 
gravitation was discussed in Refs. mil 

The bilayer graphene in a simplest approximation can 
be considered as a zero-gap semiconductor with parabolic 
touching of the electron and hole bands described by the 
single-particle Hamiltonia n 1 8 > 1 9 



H = 





(p x + iPyf /2m 



{Px - iPy) /2m 




(5) 



where pi — —ihd/dxi — Ai are electron momenta 
operators and m ~ 0.054m e is the effective mass, m e 
being the free-electron mass. This description is accurate 
at the energy scale larger than few meV, otherwise 
a more complicated picture including trigonal warping 
takes placed; we will restrict ourselves only by the case of 
not too small doping when the approximate Hamiltonian 
([5]) works. Two components of the wave function are 
originated from crystallographic structure of graphite 
sheets with two carbon atoms in the sheet per elementary 
cell. There are two touching points per Brillouin zone, 
K and K' . For smooth enough external potential, no 
Umklapp processes between these points are allowed and 
thus they can be considered independently. 

We will proof that the zero-energy states in the 
case of bilayer found by exact solution for the case of 
homogeneous magnetic field^ are topologically protected 
and their number is determined only by the total flux per 
sample, irrespective to whether the field is homogeneous 
or not, exactly as in the case of the single-laye r 2 ' 11 ' 13 . 



III. RESULTS AND DISCUSSION 

The proof is based on the theory of elliptic operators 
and on the Atiyah-Singer index theorem. Let us remind 
first some facts about it. 



Let A be a smooth compact manifold, E and 
E' smooth complex bundles over X (we shall use 
everywhere the word "smooth" in the sense of "infinitely 
differentiable") . Let D be a smooth linear differential 
operator of order m acting from C°°(E) to C CC {E I ) where 
C CC (E) is the vector space of smooth sections of E. Here 
smoothness of operator is regarded as smoothness of its 
coefficients in any smooth local coordinates. 

In local coordinates (x l ) on X the highest-order 
terms of D have a form J2 ° 41 '"' m ( x ) g§n ■ ■ ■ g J m ■ Let 
us consider the expression ° n ( x )£ii ■ ■ ■ £i m > £ e 
T*X, T*X being the cotangent bundle of X. It is 
independent on the choice of local coordinates and defines 
the homomorphism of vector bundles tt*E — > 7T*-E", which 
is homogeneous of degree m by £. Here n*E, ir^E' are the 
liftings of the bundles E, E' to T*X (see the commutative 
diagrams below; vertices on the diagrams are smooth 
manifolds and arrows are smooth maps). 



n*E ■ 



■E 



tt*E ■ 



a(D) 



■ir*E' 



T*X ■ 



X 




T*X 



This homomorphism ir*E — > ir*E' is called symbol cr(D) 
of differential operator D. The latter is called elliptic if 
cr(-D) is invertible outside zero section of T*X (that is 
invertible at £ 7^ in local coordinates). 

Elliptic operators have a good behavior—: if 
D: C°°(E) -> C°°(E') is a smooth elliptic operator then 

• All distributional solutions of D are smooth. 

• Ker D (the space of solutions of the equation Dip = 
0) and CokerD (the factor-space of C°°(E') by the 
image {Dip} of D) are finite dimensional. 

• index D — dim Ker D — dim Coker D depends only 
on the symbol of D; moreover, index D depends 
only on the homotopy class of symbol in the space 
of continuous invertible symbols of a given order. 

The rough idea of the proof is the following. We 
consider the two-periodic case (a justification of this 
choice will be discussed below). So a wave function is the 
section of complex linear bundle E over two-dimensional 
torus X. The vector potential A defines the connection 
Vj = Tj|- — jrAj on E; its flux h^ 1 J x dA defines the 
bundle E up to isomorphism. 

The differential operators p x ±ip v are elliptic operators 
acting on C°°{E). They are conjugated so the co-kernel 
of p x + ip y is isomorphic to the kernel of p x — ip y and 
vice versa. Therefore the difference of dimensions of the 
kernels of p x + ip y and p x — ip y is equal to the index of 
operator p x + ip v . The same is valid for the squares of 
these operators thus dimKer^ + ip y ) 2 — dimKer^ — 
ipy) 2 = index(p x +ip y ) 2 . 

The desired result follows just from the fact that the 
index of the composition of elliptic operators is equal to 
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the sum of their indices 3 - 3 -, so we have index^ + ip y ) 2 = 
2 index (p x + ip y ). 

According to the Atiyah-Singer theorem the index of 
operator p x + ip y depends on only the symbol of p x + ip y 
determined by the integer number N = (2ttK)~ 1 j x dA 
and does not depend on the choice of the field A for 
a given number of this integral. We obtain index(p x + 
ip y ) — N from purely topological considerations, 
replacing operator p x + ip y by other operator with the 
same symbol and known index. In physical terms, N is 
the total flux of the (pseudo) magnetic field per torus in 
the units of the flux quantum. 

It is worth to stress that the "vector potentials" A are 
assumed to be in our proof not very smooth but just 
continuously differentiable which makes the result rather 
general. 

The choice of the torus can be justified by standard 
arguments used at the introduction of the Born - von 
Karman periodic boundary conditions in solid state 
theory 34 . Namely, the total number of zero-energy modes, 
assuming that N ^ is proportional to the total number 
of atoms in the sample, Nq. At the same time, if one 
replaces "realistic" boundary conditions by the periodic 
ones the total density of states can be changed by a 
quantity proportional to the number of edge atoms, that 
is, t/No. This means that the total number of states 
with the energy close to zero should be, in the limit of 
large crystallite, correctly described by the periodic case, 
that is, the case of torus. Note that the torus has zero 
Gaussian curvature which physically means absence of 
topological defects, such as disclinations (pentagons of 
heptagons in the original hexagonal lattice) 31 ' 32 . 

The physical consequences are straightforward: as well 
as for the case of the single-layer graphene, for the case 
of bilayer (i) corrugations can result in the appearance of 
the mid-gap states^ and (ii) pseudomagnetic fields due 
to corrugations will not broaden the zero-energy Landau 
level in the case of quantum Hall effected. It would be very 
interesting to check experimentally the second statement 
by measurements of the quantum Hall activation gaps 
for the bilayer graphene, similar to Ref. |3(J for the case 
of single layer. 

A mathematical proof of the statement is presented in 
the Appendix. 
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Appendix 

Let r be a lattice in the two-dimensional Euclidean 
space R 2 , X = R 2 mod T be the two-dimensional torus, 
£ be a smooth linear complex vector bundle over X 
with structure group 17(1). Let we have C 1 -connection 
on E, that is corresponding covariant derivatives are 
written as V x — ^ - iA x , V y — ^ - iA y , where A x , 
A y are continuously differentiable real functions in local 
coordinates (x, y) on X (these coordinates we choose as 
usual coordinates on universal covering K 2 of A, so they 
are defined up to addition of vectors from T) . 

Although the connection form A = A x dx + A y dy 
depends on the choice of local coordinates, its curvature 
dA is globally defined on X. We can integrate dA over 
X; this integral depends only on isomorphism class of E 
and does not depend on the choice of connection on E. 
Let 

N{E) = (27T)- 1 f dA; (6) 
Jx 

this number must be integer. 

Following Ref. [TH, we consider Hilbert spaces H S (E) 
of those distributional sections u of E for which Du e 
L 2 (X) for all differential operators D: C°°(P) -> 
C°°(lx) with smooth coefficients, and of order < s. 
Here C°°(E) is the space of smooth sections of E, lx is 
the trivial linear complex bundle over X. The Hermitian 
product in H S (E) can be defined at s = as (u,v) = 
J x (it, v) dxdy, at s > - as (u, v) s = J x (A s u, v) dxdy. 
Here A = 1 + D*D, D : C°°(E) -> C°°{E <g> T*X) 
is the covariant derivative given by some fixed smooth 
connection on E (precise choose of this connection is 
irrelevant for our aims). 

We can consider differential operators P ± = V x ± iV y 
as continuous linear operators from H S (E) to H S -\{E) 
at s < 2; let us denote these linear operators as P t ,. 
Similarly, we can consider differential operators Q ± = 
(Vz ± iV y ) 2 as continuous linear operators Qf from 
H 2 (E) to H {E). 

Theorem. 

dimKerP+ - dimKerPf = N(E), 

1 (7) 
dimKerQ+ - dimKerQ^ = 2N(E). 

Remark. If the connection (that is the functions 
Ax, A y ) is smooth then all distributional solutions of 
the operators P^, Qf are also smooth^, and for the 
differential operators , acting on C°° (E) we have 
from Eq.(0 

dimKcrP + - dimKerP" = N(E), 
dimKcrQ+ - dimKerQ" = 2N{E). 

Proof. Note that in our case when X is the two- 
dimensional torus the cotangent bundle it: T*X — > X 
is trivial two-dimensional real bundle over X and can 
be identified with the trivial linear complex bundle 
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X x C -> X. So lift of the bundle £? over X to the 
bundle ir*E over T*A can be identified with linear 
complex bundle ExC^XxC. At this identification 
symbols cr ± : it*E — > 71% P of the operators P ± become 
the following form: for e € £7, £ € C we have <r + (e,£) = 
(? e )0; ^"(ejC) = (£ e j£)> that is the fiber over a point of 
T*X is multiplied by the complex number corresponding 
to this cotangent vector in the case <r + , and on the 
conjugate to this complex number in the case o~ . 

The composition a + a~ : tt^E — > tt*E, cr + er~(e, £) — 
(l£| 2e j£) coincides with the identity id: n*E — > n*E on 
the unit sphere bundle of T*X. So [er + ] + [<r~] = 0, 
where [a] is the class of a in the group K (T*X) where K 
denotes if -theory with compact supports (a description 
of this variant of if -theory is contained in Ref . [T3) . 

Applying "topological index", that is, homomorphism 
index: K(T*X) — > Z constructed by Atiyah and 
Singer—, to this equality, we get 

index[er + ] + index[cr~] = 0. 

The operators P±: H S (E) -> if s _i(.E) are Fredholm 
since symbols a ± are invertible outside the zero section of 
T*X (Ref.. 3). index = dim Ker - dim Coker P± 
depends only on [cr ± ] and are independent of the choice 
of s < 2 and connection field A (but of course they 
depend on N(E), which define the isomorphism class of 
E): index P± = index^] (Ref. El). 

Note that at s > 1 for u,v £ H S (E) we have 
uv £ H s (lx). So J x (ww) K da;dy = f x (uv) y dxdy = 0, 
and (P 1 + u,w) + (tt, Pfu) = for any u, v £ H\(E). 
Identifying Coker Pj^ with the orthogonal complement 
of ImP^ in Hq(E), we obtain 

Ker Pf = (Coker Pf) n H X (E) 
Ker Pf = (Coker Pf) n ffi(-E) 
Hence 

index Pf < dim Ker Pf — dim Ker Pf 
index Pf" < dim Ker Pf" — dim Ker Pf 1 " 

where every of these inequalities become equality if 
the co-kernel of the corresponding operator contains in 
Hi(E). However, index Pf + index Pf" = index[er + ] + 
index[cr~] = 0, so both inequalities should be equalities, 
and we obtain Coker P^ C Hi(E), and KerP* = 
Coker Pf. Hence, 

dim Ker Pf*" — dim Ker Pf = index Pf 1 " = index[cr + ]. (8) 

Repeating this consideration almost literally for the 
operators Qf — (W x ±iW y ) 2 — P^P^ 1 acting from 
Hz(E) to Ho(E), and using the fact that index of the 
composition of Fredholm operators is equal to the sum 
of their indices^, we have 

Coker Q± C H 2 (E), 

KerQf ^ Coker Qf, 

dim Ker — dim Ker Q^" = index Q\ = index(P 1 + P 2 + ) = 

= index Pf + index P 2 + = 2 index[cr + ] . (9) 



We present below an explicit calculation of the value of 
index[cr + ] based on one famous theorem from algebraic 
geometry. But for more clarity we start with a simple 
reasoning showing the proportionality of index[<7 + ] to 
N{E). 

Let us see on the construction of complex bundles 
over X. Let F be a £/(n)-vector bundle over X. Cut 
out the disk B 2 from the torus X. Since the disk is 
contractible, the restriction of F to B 2 is trivial. X — B 2 
is homotopically equivalent to the wedge product of two 
circles, and U(n) is connected, so the restriction of F 
to X — B 2 is trivial, too. Thus, the isomorphism class 
of the bundle F is uniquely defined by its dimension 
n and by the homotopical class of the map ip: S 1 — 
dB 2 — > U(n) gluing together two these trivial bundles. 
This homotopical class is defined by the degree of the 
map det-y: S 1 — > U(l), where det: U(n) — > U(l) is 
the determinant homomorphism. The sum of degrees 
corresponds to the Whitney sum of a bundles over X, so 
K{X) = Z © 1. Here K(X) is Abelian group generated 
by elements [F] with relations [F F'} = [F] + [F'] for 
all complex bundles F, F' over X; detailed description 
of if -theory is contained in Refl35l. 

Particularly, the isomorphism class of a linear bundle 
E is defined by the integer N(E) — deg(det -tp) = 
(2TT)- 1 J x dA, and for the class [E] of E in K{X) we 
have 

[E] — l = N(E)([Ei] — 1), (10) 

where E\ is a linear bundle over X for which N(Ei) = 1. 

If E is trivial then choosing trivial connection we 
obtain that Ker P + is the space of holomorphic functions 
on torus and KerP - is the space of anti-holomorphic 
functions on torus. Both these spaces contain only 
constants and are 1-dimensional, so index P + = in this 
case. Taking into account that a + is the image of [E] € 
K(X) at the Thorn isomorphism K(X) -> K(X x C) 
(the description of this isomorphism see in Ref. , from 
(fTO)! we obtain that index[er + ] is proportional to N(E). 

To calculate the coefficient of this proportionality, 
moreover, to calculate the value of index[cr + ], let us 
replace E by the other bundle of the same class in K(X), 
and replace Pf by the other operator of the same symbol 
class in K(T*X) (the index of the operator does not 
change at such a replacement). 

Consider the torus X as algebraic curve, with local 
complex coordinate z — x + iy. Choose holomorphic line 
bundle F over X isomorphic to E in smooth category, 
that is such that N(F) = N(E) (for example, we 
can take divisor on X consisting of a point zq € X 
of the multiplicity N(E), and turn from the divisor 
to corresponding holomorphic line bundle by the way 
described in Ref. l36h . 

Consider now the differential operator 8: C 00 ^ ' ) — > 
C°°(J*.i), 8 = £ = £ + ifa where T^ k is the 
bundle of differential forms on X of type (0, k) with 
coefficients in F. Since the complex cotangent bundle 
of X is trivial and linear, we have J- 0,0 = J- 0,1 = 
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F, so the symbol of the operator d is coincide with 
a + . Hence index[cr + ] = index 8. However, index d is 
equal to the Euler characteristic x(-X", F) of the sheaf 
of germs of holomorphic sections of F. We can compute 
X(X, F) using the Riemann-Roch-Hirzebruch theorem^. 
For a curve X and linear bundle F this theorem yield 
X{X, F) = cx(F)[X] + 1 — g, where g is the genius of the 
curve X and c±(F) € H 2 (X; Z) is the first Chern class of 



the bundle F. In our case g = l, d{F)[X] = q(E)[X] = 
(27r) _1 J x dA — N(E), so we obtain the final formula 

index[cr+] = N(E). (11) 

Substituting ifTTj) to (fij-fljlj, we obtain the assertion of 
the Theorem. 
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